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Abstract 

This paper presents the shape optimization of a flat-type arborescent fluid distributor for the pur- 
pose of process intensification. A shape optimization algorithm based on the lattice Boltzmann 
method (LBM) is proposed with the objective of decreasing the flow resistance of such distribu- 
tor at the constraint of constant fluid volume. Prototypes of the initial distributor as well as the 
optimized one are designed. Fluid distribution and hydraulic characteristics of these distributors 
are investigated numerically. Results show that the pressure drop of the optimized distributor 
is between 15.9% and 25.1% lower than that of the initial reference while keeping a uniform 
flow distribution, demonstrating the process intensification in fluid distributor, and suggesting 
the interests of the proposed optimization algorithm in engineering optimal design. 

Keywords: Lattice Boltzmann method. Fluid distributor. Shape optimization. Heuristic 
optimality criterion. Flow distribution. Pressure drop 



1. Introduction 

Delivering and distributing flows of one or different fluids precisely onto a given surface or 
into a given volume is an important issue in many unit operations of process engineering. For 
example, homogeneous flow distribution is generally required for multi-channel heat-exchangers 
[1-3]; for solar collectors [4-6]; for trickle bed, catalytic packed bed or multi-channel reactors 
[7-12]; for catalytic monoliths [13,14] or for bubble column [15,16]. 

It has been well recognized that the flow distribution uniformity is generally associated with 
proper design of entrance and exit manifolds, or the configuration of fluid distributors/collectors 
[17,18]. Many researches have focused on the design and optimization of fluid distributors/collectors, 
in particular, arborescent architectures have been given much attention; whether from one point 
to a line [9,12,19,20] or to the periphery of a disc [21-23]; to a surface [24-27] or to a vol- 
ume [28,29]. The basic architectures of multi-scale channel networks build on analogies with 
living organisms (e.g. the lung or the vascular system). Design and optimization procedures 
have been proposed by Tondeur and his coworkers [26,27,29,30] to minimize the pressure drop 
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(viscous dissipation) under constraints of uniform irrigation and constant void volume based on 
the principle of "equipartition properties", which is in close connection to the constructal theory 
developed by Bejan [31-33]. Analytical, numerical and experimental investigations carried out 
by Luo and her coworkers [29,34,35] imply that arborescent structures that may guarantee the 
flow distribution uniformity between the channels will induce higher pressure drops, mainly be- 
cause of the so-called "minor losses" due to numerous singularities of the complex multi-scale 
structure (bifurcation, elbow, etc.)- These losses are actually not so minor As a result, for the 
design of an arborescent fluid distributor, can the analytical optimization be further improved, 
particularly for the profiling of the junctions in complex structures? 

In the last decades, numerical techniques have been developed to solve the problem of flow 
shape optimization. Earlier work has studied various aspects of shape/topology optimization of 
fluid flow [36-38]. Noteworthy is the work of Errera and Bejan [39] which proved that the den- 
dritic patterns formed by low-resistance channels in a river drainage basin can be deduced from 
the constrained minimization of global resistance in area-to-point flow. Moos and his coworkers 
[40,41] also proposed a procedure for topology optimization of fluid flow based on the principle 
that the fluid flow always searches for the best way under given constraints in a predefined space 
by itself. 

So emerges the idea that classic optimization approaches could benefit of numerical com- 
puting to overcome the analytical difficulty of expressing objective function and constraints in 
complex geometries. In that context, we have developed a lattice Boltzmann method (LBM) 
shape optimization algorithm for minimizing the resistance of flow structures at the constraint 
of fixed void volume for fluid flow [42]. Numerical examples of a right angle elbow and a T- 
junction show that this algorithm can optimize the flow shape with significantly reduced flow 
resistance. In the present study, this algorithm is applied to optimize the structure of a flat-type 
arborescent distributor investigated in reference [27]. Numerical experiments are also carried 
out in parallel for comparing flow distribution and hydrodynamic characteristics of the distribu- 
tors. Results obtained are useful for further validation of our optimization algorithm, showing its 
promising application in processes dealing with fluid. 

2. Numerical algorithm 

LBM, originally proposed by McNamara and Zantti (1988) [43] as a smoothed alternative to 
lattice gas automata (LGA), is an efficient second-order Navier-Stokes solver capable of solving 
various systems for hydrodynamics owing to its algorithmic simplicity, explicit solution of par- 
ticle distribution functions, natural parallelism and easy boundary treatment [44]. The objective 
of this algorithm is to optimize the shape of fluid flow by minimizing the flow resistance (thus 
the pressure drop) under the constraint of constant void volume for fluid flow. To do that, a 2D 
simulation domain including both fluid phase and solid phase is uniformly divided into elemental 
square cells. LBM is utilized as an underlying Navier-Stokes solver to calculate the flow flied. 
Compared to the traditional CFD methods based on discretizations of macroscopic continuum 
equations, the kinetic nature and local dynamics of LBM make it more adaptive to complex 
boundaries and parallel computing [44]. In addition, the theoretical basis ("ceU" expression) of 
the LBM corresponds very well to the downstream cells' position exchanging procedure. That is 
the reason why we choose the LBM as the "pre treatment". 

To simplify the numerical algorithm, following assumptions are made: 

- Steady flow pattern; No-slip condition at the wall. 
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Figure 1: A two-dimensional (2D), 9-velocity D2Q9 lattice model 

- 2-D simulation, i.e. infinite channel depth; only consider the friction between the fluid and 
the solid walls at both sides; negligible gravity effect; 

- Isothermal operating condition; isotropic and homogeneous physical properties of solid 
materials; constant physical properties of working fluid. 

2.1. Standard lattice Boltzmann formulations for fluid flows 

In standard lattice Boltzmann method, the macroscopic equations of traditional CFD (the 
Navier-Stokes equations) are not solved directly, but rather the Boltzmann equation is solved on a 
discrete lattice. For the solution of incompressible fluid flow, we use the classical D2Q9 scheme 
of LBM, with three speeds (0, 1 and V2) and nine velocities (co ~ Cg) on a two-dimensional 
square lattice (Fig. 1), the evolution of the distribution function f (x, t) obeys the lattice Boltzman 
equation (LBE) [44]: 

f (X + c,Af, t + At)- f (X, f) = -i [f - f^") , /=0,1,2 ■ ■ ■ 8 (1) 

Herec,= {cos [(/ - l)7r/2] , sin [(/ - l)7r/2]) for/= 1 ~ 4, c,= V2 {cos [(/ - 5) ;r/2] , sin [(/ - 5)7r/2]) 
for 5 ~ 8 and Co=0, and r = 3v -t- y is the relaxation time related to kinematic viscosity v and 
the discrete time step At, and f^'* is the equilibrium distribution function defined as [45], 

^q. ^ c,-u l(cru)^ lu-u\ . 

where the weights are given by (jJo= 4/9, (^, = 1/9 for /= 1 ~ 4, (^, = 1/36 for /= 5 ~ 8, and Cj is the 
speed of sound and equals to V3/3 in lattice unit, p and u are density and velocity, respectively. 

The macroscopic variables such as fluid density and momentum are related to the distribution 
function and obtained as: 

p (X, t)^Y^f (x, f), pu (x, = 2 Cif (x, f) (3) 

/ / 

Through the Chapman-Enskog procedure, the incompressible Navier-Stokes equations can 
be obtained from the LBE in the limit of small Mach number [46-48], as the following: 

|+V.(pu) = 

(4) 

^ -H V . (puu) = -V(pc2) + V . (pvVu) -hF 




where the pressure is given by: 



(5) 



2.2. Mutual replacements between fluid and solid cells 

We pay close attention to the dynamic interaction at fluid-solid interface and propose a heuris- 
tic optimaUty criterion, i.e. the viscous stress for solid cells: 

and dynamic pressure for fluid cells: 

q = \p\\nf (7) 

where Ax and Ay are the lattice spacing (Ax = Ay in the LBM), u and v are velocity in x and y 
direction, respectively. 

In our optimized algorithm, solid cells at fluid-sohd interface suffering from larger viscous 
stress from its neighboring fluid cells will vanish and be replaced by fluid ceUs. Likewise, fluid 
cells at the fluid-solid interface having lower dynamic pressure are replaced by solid cells, in 
order to ehminate the "dead zones" in the fluid domain so as to effectively make use of a fixed 
void volume for fluid flow. Note that the equal number of fluid cells and solid ceUs is targeted to 
balance the void volume occupied by fluid. The mutual replacement of fluid and solid cells will 
create a new flow shape that differs from the initial one. This cells' position exchange process 
in our algorithm can be considered as a "generalized Cellular Automaton" and used the term 
"CA" for short. Then the flow field of this new shape will again be calculated by LBM for 
the recurrence by a CA procedure. Step by step, the shape of fluid flow evolves towards the 
final shape with homogenized dynamic pressure at fluid-sohd interface and reduced total flow 
resistance. 

From the view point of fluvial dynamics, our algorithm can be thought of as a mimicry of 
natural behavior in river channels where the surface is eroded at the points of maximum shear 
stress and the sand is deposited at the points of minimum dynamic pressure [42]. 

2.3. Numerical implementations 

The optimization procedure is described in detail using the flow chart shown in Fig. 2. 

(1) Input the initial data such as the size and initial shape of the simulation domain (solid 
phase, fluid phase), the specified boundary conditions (fluid nature, velocity profile, pres- 
sure, etc.). 

(2) An exact flow field is calculated by LBM. 

(3) At the fluid-solid interface, a number of fluid cells having the lowest dynamic pressure and 
the same number of solid cells suffering the largest viscous stress will be identified, and 
their positions will be exchanged, thus creating a new shape. 

(4) Reinitialize the basic properties of fluid flow and boundary for the new shape. Recalculate 
the exact flow field by LBM. 

(5) Check the stable tolerance of the algorithm. If the tolerance is satisfied, then the heuristic 
procedure is terminated, and the results are exported. Otherwise, the procedure goes back 
to Step 3 for recurrence. The result is considered to be stable when the pressure drop Ap 
across the system tends to extremum. 
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Figure 2: Flow chart of LBM based shape optimization algorithm 
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Figure 3: A simple numerical example of LBM based shape optimization algorithm 



We have tested an elemental T junction using our optimization algorithm, as shown in Fig. 
3. It can be observed that T junction evolves into a Y shape junction with about 75% decrease 
of A/7. It is expected that the total pressure drop of the arborescent structure having successive T 
junctions could also be further reduced by applying our algorithm. 

3. Design and optimization of flat type distributor 

3.1. Conventional flat type distributor 

Conventional flat type distributor comprises a simple cuboid distributor body having a fluid 
port at one face and a distribution surface at the opposite face, with the distribution surface 
having a plurality of uniformly spaced distribution openings. We choose the conventional flat 
type distributor as our comparison target. 

3.2. Arborescent distributor 

A flat type arborescent distributor is also designed for comparison. As shown in Fig. 4, the 
tree-shaped structure is induced by four generations of T-bifurcation or division, the number of 
outlets being 2"* - 16, uniformly distributed in a square surface. The four generations or scales 
are indexed from to 4, with for the single inlet channel and 4 for the smallest channels. Each 
outlet is 20 mm apart from its horizontal and vertical neighbors. For this reason, the channel 
length at each scale is employed as the "local constraint", which obeys: 



The cross-section of the distributor's channels is rectangle. The optimal scale ratio for the 
channel width at each scale is obtained by an optimization that accounts for both viscous 
dissipation and total pore volume, with the assumption of established laminar Poiseuille flow 
and neglecting the effects of flow singularities [26]. Detail optimization procedure may be found 
in the reference [27], and the scaling law for channel width: 



h-h- 20 mm; - U - 10 mm; 



(8) 




(h + Wk+i)(3h + Wk+i) 4 
6 



(9) 




Figure 4: The arborescent structure of the flat-type distributor 



Table 1: Dimensions of the arborescent distributor 



Channel 



Index 


Number n 


Width w 


Depth h 


Length / 


1 


2 


2.00 mm 


10 mm 


20 mm 


2 


4 


2.91 mm 


10 mm 


20 mm 


3 


8 


4.27 mm 


10 mm 


10 mm 


4 


16 


6.35 mm 


1 mm 


1 mm 



Where h is the depth of the channels, which is considered as constant (10 mm) in a flat-type 
distributor. In this work, W4 is set to be 2 mm for the fabrication tolerance to be acceptable. 
With this specification, the dimensions for the bifurcated channels can be calculated using Eq. 
(9) and are listed in Table 1 . Good flow distribution by this type of arborescent distributor has 
been reported because of its path symmetry that provides equivalent hydrauhcs characteristics, 
i.e. equivalent flow rate, equivalent time of passage and equivalent pressure drop [27,30]. 

3.3. Optimization by LBM based algorithm 

Based on the arborescent structure shown in Fig. 4, the entire computation domain is uni- 
formly divided into elemental square sub-domains, which are categorized into fluid cells, solid 
cells, inlet cells and outlet cells (Fig. 5). Considering the symmetry of simulation domain, half 
of the simulation domain is computed. The simulation domain is divided into 180 x 360 cells. 
The width of inlet and that of outlet are Win = 11 cells and Wout = 52 cells, respectively. The 
parabolic velocity profiles are imposed at inlets where the flow rates are identical and the Neu- 
mann boundary condition [49] at the outlet is enforced. The waU boundary condition given in 
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Figure 5: Cellulai' expression and boundary condition 



this paper is the bounce -back scheme, which is adaptable, robust and easy to be implemented 
for fluid flow in complicated geometries. The fluid density and kinematic viscosity are po- 1.0 
and V- 0.04, respectively. The initial condition is the equilibrium distribution of the inputs over 
the whole computational domain, using a constant density po= 1-0 and zero velocity. The cross- 
sectional average fluid velocities of inlets are Ui„ = 0.08, corresponding to a Reynolds number 
Re^^ = ^ 110.08/0.04 = 22. 

The evolution of shape and flow field at different optimization time steps is shown in Fig. 6. 
With the algorithm proceeds, the flow singularity (right angle) gradually disappears and becomes 
smoother, implying the alleviation of singularity effect. Finally, it reaches a relative steady state, 
i.e. the shape of ffuid flow stays almost unchanged with increasing time step. It can be observed 
that the streamline of velocity profile is weU-kept at the final shape so that the flow turns slowly 
and continuously rather than abruptly. 

It should be pointed out that what we optimized actually is an arborescent fluid collector 
for low Reynolds number applications. Standard LBM is developed for laminar flows at low 
Reynolds numbers, not applicable for problems with small kinematic viscosity as they associ- 
ated with high Reynolds numbers, or being turbulent in many practical applications. Therefore, 
the incorporation of turbulence model into LBM together with heuristic optimality strategy to 
address practical systems at high Reynolds numbers will be well considered in our future work. 
Generally, there are two alternative ways to model turbulence: introducing the Reynolds equa- 
tion and turbulent stress as in - e model or using space-filtered governing equation and large 
eddy simulation (LES) with subgrid-scale stress model for the unresolved scale stress. 

3.4. Prototypes set-up 

Figure 7 shows the inner structures of conventional distributor, arborescent distributor and 
optimized distributor, respectively. The top has a thickness of 5 mm with a hole of 10 mm in the 
center connected to the inlet tube 10 mm). The middle one having different structures . The 
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Pore structure 



(a) Conventional distributor 




(b) Aborescent distributor (c) Optimized distributor 



Figure 7: Pore stmcture of the flat-type distributors 

bottom plate is 5 mm thick and has 16 holes ((iout=2 mm) located at the 16 ends of the intemal 
channels on the middle plate. Connecting to the holes are 16 tubes having an internal diameter 
of 2 mm and a length of 10 mm. They serve as the outlets of the distributor 

4. Numerical simulations 

The distributors together with the initial one, optimized by constructal laws and by LBM 
optimality algorithm are assigned respectively as conventional distributor, arborescent distributor 
and optimized distributor, have been tested numerically, in order to investigate and compare their 
flow distribution and hydraulic performances. 

4.1. CFD simulation 

For the CFD simulation, the following assumptions were made: 

- Fluid flow is incompressible, isothermal, Newtonian, in steady-state; the effect of viscous 
heating is negligible. 

- The physical properties of the fluid (liquid water at 20°C) are: density p - 998.23 kg ■ m"^; 
viscosity ji - 0.00101 kg ■ m"' ■ s"'. The operating pressure is set at 101325 Pa. 

The geometry model was built up using Hex-wedge elements by software GAMBIT 2.3.16. 
Note that with the symmetric assumption, one quarter of the real object was adopted as the model 
for the purposes of lessening the computational burden. Computational grids with 5.0x10^ nodes 
were selectively refined in some local place where parametric variation was severe. 
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CFD commercial software FLUENT 6.3.16, of ANSYS, Inc., USA, was used here. Navier- 
Stokes equations were solved in 3D by standard k - e segregated turbulent solver with stan- 
dard wall functions for near-wall treatment.The SIMPLE algorithm was employed for pressure- 
velocity coupling. The momentum equations were solved with second-order upwind differencing 
while the pressure terms were discretized using the standard scheme. At the inlet port, a constant 
velocity profile normal to the entry surface was used as initial boundary condition. The outlet 
ports were defined as pressure-outlet, with a zero relative pressure. A no-slip condition was 
imposed at the walls. The calculations were performed in double precision. The solution was 
considered to be converged when (i) the mass flow rate at each channel and the inlet static pres- 
sure were constant from one iteration to the next (less than 0.1% variation) and (ii) the sums of the 
normalized residuals for control equations were all within the order of magnitude of 1 .0 x 10"^\ 

The simulations were conducted under the inlet water velocity ranging from 0.01 to 3.0 
m ■ s"', corresponding to a Reynolds number within the range of 99.4~29821.0in the inlet chan- 
nel and an average Reynolds number within the range of 140.4~46679.0 in the outlet channel. 

4.2. Measurements of flow distribution uniformity 

Two dimensionless parameters, namely the maximal flow-rate ratio 6 and the maldistribution 
factor Dg, were used to evaluate the flow uniformity of distributors: 

e = ^ (10) 



1 ^ 
1 , m, 



1 j-t »Jave 



Where, A' stands for the sampling channels' number and Wave the mean flow rate of channels: 



1 ^ 

Wave = X (l^) 
1=1 



5. Discussion and Conclusion 

5.7. Flow distribution 

CFD simulations have been performed to investigate the flow distribution performance of the 
initial and optimized distributors. Showing in Fig. 8 are the variation of maximal flow-rate ratio 
and the maldistribution factor for three distributors as a function of average Reynolds number at 
inlet. It can be observed that under the operating conditions examined, influence of flow-rate on 
the flow distribution uniformity for conventional distributor is large while this influence is very 
small for both arborescent distributor and optimized distributor. Both the maximal flow-rate ratio 
and maldistribution factor for conventional distributor are larger than those of both arborescent 
distributor and optimized distributor and decrease with the increasing average Reynolds numbers 
at inlet, implying that the fluid distribution performance of both arborescent distributor and op- 
timized distributor are far better than that of conventional distributor. The flow maldistribution 
in both arborescent distributor and optimized distributor is within a very low level: the maximal 
flow-rate ratio less than 1.05 and the maldistribution factor less than 0.015. Note that the flow 
uniformity is a bit better for the arborescent distributor than the optimized one. The influence 
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of the optimized algorithm on fluid distribution comes from two factors. On one hand, the op- 
timized channel shape is not as symmetrical as before, at a bifurcation point, more fluid will go 
to the sub-branch with a small turning angle owing to the inertia of the fluid. And on the other 
hand the resolution of the square lattice for the LBM based optimized algorithm may be not high 
enough. The relative coarse boundary generated for the optimized construct would be transposed 
as input for geometrical model set-up for simulation, which may cause slight geometry differ- 
ences and consequently very slight increase of flow maldistribution (ADg < 1%). Nevertheless, 
we emphasize that both distributors can guarantee excellent flow distribution uniformity under 
our operating conditions, with a large improvement with respect to the conventional case. 

5.2. Pressure drop 

Figure 9 shows the contour of static pressure for the reference construct and for the optimized 
one. It can be clearly observed that there is a large pressure drop at the outlet of the distributor, 
which is more than 50% of the total pressure drop. Besides the vertical inlet channel (scale 
0), the pressure drops in the three distributors are generally located at the distributing network, 
and the frictional pressure drop in outlet parallel tubes is negligible. Examining the arborescent 
distributor, we observe that pressure drops are generally generated in T-type junctions where 
flow spUts, because of strong viscous friction induced by the crash between the fluids and wall 
located in front of the flowing direction. This phenomenon has already observed and reported in 
our earlier work [34, 35]. Contrarily, the eflFects of singularities in the optimized case are largely 
eliminated because of the smooth transitions instead of the right angles in the Y-type junctions. 

Figure 10 shows the comparison of pressure drop for the three distributors as a function of 
average Reynolds number at inlet. It is should be pointed out that the pressure drop of distributor 
is referred as the pressure difl'erence between inlet tube connecting to the hole at top of distributor 
and 16 small outlet tubes connecting to the holes at bottom of distributor The inlet static pressure 
of each distributor is within an acceptable range (less than 0.35 bar) under our tested conditions. 
The pressure drop of each distributor rapidly increases with increasing flow-rate. At the same 
flow-rate, the pressure drop of conventional distributor is far smaller than that of both arborescent 
distributor and optimized distributor. It should be noted that the pressure drop of arborescent 
distributor is consistently higher than that of optimized one because of its relatively higher flow 
resistance. The pressure drop in the distributor is reduced from 15.9% to 25.1% when average 
Reynolds number at the inlet increases from 99.4 to 29821.0, implying the good performance of 
our LBM based shape optimization algorithm in design and optimization of fluidic devises. 

6. Concluding remarks 

A flat-type arborescent distributor is optimized using a LBM based shape optimization al- 
gorithm developed by us. Distributor prototypes with conventional, arborescent and optimized 
arborescent structures are designed and investigated numerically. Numerical results show that 
both optimized distributor and arborescent distributor can guarantee excellence flow distribu- 
tion uniformity while the conventional distributor increases maldistribution, the pressure drop of 
the conventional distributor is far lower than that of both arborescent distributor and optimized 
distributor and the pressure drop of the optimized flat-type distributor can be reduced between 
about 15.9% and 25.1% while keeping a uniform flow distribution, demonstrating the good per- 
formance of our optimization algorithm in design and optimization of complex structures of fluid 
flow. 
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Figure 8: Flow distribution uniformity of the three distributors as a function of average Reynolds number at inlet, (a) 
maximal flow-rate ratio; (b) maldistribution factor 



13 



I 



4,996+03 ^ 
4.0Se+03 

3.12^.03 . \ 

125e+03 ^^PVW ^ 

3.10e-f02 ^^HF ^i^k 

-6.25e+02 41^^ X. ^ 

-1.56e+03 ^ 
-2.50a+03 

-3,43e+03 ^ . ,.Y 
-4,37e+03 ^ 



(a) Conventional distributor 




I 



.56e+D3 
50e+03 
.43e+03 ; 
.37e+03 
,30e+03 




(b) Arborescent distributor 



(c) Optimized distributor 



Figure 9: Contour of static pressure for three distributors (Inlet velocity: 2.0 m • s ; average Re at inlet: 19881.0) 
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Figure 10: Pressure drop versus average Re at inlet for three distributors (inlet velocity range : 0.01 -3.0m-s-') 
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Two different numerical methods, i.e. the discrete-velocity Boltzmann equation based CFD 
method (LBM) and the traditional Navier-Stokes based CFD method (Fluent) are used in this 
study. Indeed, our shape optimization algorithm is proposed as a method of searching for op- 
timal architecture with minimal constraints. The limitation comes essentially from the initial 
morphology: there is no "freedom" to "morph" in the traditional Navier-Stokes based CFD as 
the geometry is fixed; there is a Umited "freedom" to "morph" in classic analytical approaches 
(the establishment of scaUng relations); and there is a relatively large "freedom" to "morph" with 
LBM approaches. 

In general, the optimized distributor could find a balance between arborescent distributor 
and conventional distributor in both uniform distribution and pressure loss, and the pressure 
drop decrease is not as significant as expected with respect to the single T shape case (Fig. 
3). Indeed, it would be more convenient by rounding the sharp comers in engineering cases. 
However, our algorithm actually tries to give an extreme of the pressure drop reduction. There 
might be two main reasons for this limited reduction. Firstly, this is due to the fact that the 
initial arborescent structure taken as the starting point for our optimization algorithm is already 
analytically optimized, as indicated by the established scaling relation of Eq. (9). What we 
searched for by applying the LBM based shape optimization algorithm is actually the "further 
improvement" and the results obtained seem reasonable and encouraging. Secondly, the limiting 
factor lies in the 2D nature of our current optimization algorithm, i.e. the height of the channels 
is assumed infinite. In fact, the friction between the flow and the upper/bottom walls of channels 
with rectangular section, and the entrance elTect (viscous dissipation in single inlet channel and in 
the first splitting T junction) contribute much to the total pressure drop in the flat-type distributor, 
as indicated by CFD simulation results (Fig. 9). A 3D optimization algorithm is expected by 
employing 3D lattice Boltzmann models, taking the influence of channel's cross-section shape 
(circular instead of rectangular) into account, for the optimization of 3D distributors presented in 
the references [26,29]. 

It should be pointed out that the optimized structure we obtained is a flat type arborescent 
fluid collector for low Reynolds number appUcations. However, CFD simulation results show 
that the performance improvement is obvious, even under turbulent flow conditions, indicating 
the robustness of the optimized architecture obtained using our algorithm. Incorporating proper 
turbulence model into LBM together with heuristic optimaUty strategy for high Reynolds number 
conditions will be another direction of our future work. Note that the decrease in total pressure 
drop might be even more evident for the use of collecting fluid flow (fluid collector), which is 
currently not able to be tested experimentally. 

Finally, a further goal is to develop our algorithm for the optimization of fluid flow structures 
to optimize the performances of heat exchangers or reactors . In that case, it's clear that other 
heuristic optimahty criteria that account for fluid flow, heat transfer or reaction kinetics should 
be introduced. 

Nomenclature 

Ap pressure drop, L"^MT"^ 
At time interval, T 
Ax lattice spacing in x direction, L 
Ay lattice spacing in y direction, L 

15 



/Wave 


average flow-rate, MT"' 


'^max 


maximum flow-rate between channels, MT~' 




minimum flow-rate between channels, MT~^ 


J" 


fluid viscosity, L~^MT"' 


V 


kinematic viscosity, L^T~' 




lattice weighting factor, dimensionless 


u 


cross-sectional average fluid velocity, LT"' 


Wout 


cross-sectional average fluid velocity of outlet, LT^' 


P 


density, L^^^M 


Po 


mean density, L^-'M 


T 


relaxation time, T 




shear stress in y direction, L~^MT"-^ 




shear stress in x direction, L~^MT"-^ 


C 


particle velocity, LT~^ 


U 


macroscopic velocity, LT"^ 


X 


particle position, L 


e 


maximal flow-rate ratio, dimensionless 


V 


velocity in y direction, LT~^ 


Cs 


speed of sound, LT"' 




maldistribution factor, dimensionless 


din 


diameter of the single inlet channel, L 


^out 


diameter of the outlet tubes, L 


fi 


particle distribution functions, dimensionless 


fT 


particle equilibrium distribution function, dimensionless 


h 


channel depth, L 


L 


channel length, L 


rik 


number of channels in scale k, dimensionless 


P 


pressure of fluid, L~^MT"^ 


1 


dynamic pressure of fluid, L'^MT"-^ 
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Re 


Reynolds number, dimensionless 


Rem 


Reynolds number at inlet, dimensionless 


t 


time step, T 


u 


velocity in x direction, LT"^ 


Min 


velocity of inlet, LT~^ 




width of inlet, L 


Wout 


width of outlet, L 


Acronyms 


CA 


cellular automata 


CFD 


computational fluid dynamics 


D2Q9 


tow-dimensional nine-velocity lattice Boltzmann model 


LBE 


lattice Boltzmann equation 


LBM 


lattice Boltzmann method 


LES 


large eddy simulation 


LGA 


Lattice Gas Automata 
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